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Abstract 

The stability of the motion of the planet satellites is considered in the model 
P_( I of the general three-body problem (Sun-planet-satellite). "Sundman surfaces" are 

constructed, by means of which the concept " Sundman stability" is formulated. The 
comparison of the Sundman stability with the results of Golubev's c^h method and 
with the Hill's classical stability in the restricted three-body problem is performed. 
The constructed Sundman stability regions in the plane of the parameters " energy - 
moment of momentum" coincide with the analogous regions obtained by Golubev's 
method, with the value {(?h)cT- 

The construction of the Sundman surfaces in the three-dimensional space of the 
specially selected coordinates xyR is carried out by means of the exact Sundman 
inequality in the general three-body problem. The determination of the singular 
points of surfaces, the regions of the possible motion and Sundman stability analysis 
are implemented. It is shown that the singular points of the Sundman surfaces in 
the coordinate space xyR lie in different planes. Sundman stability of all known 
natural satellites of planets is investigated. It is shown that a number of the natural 
satellites, that are stable according to Hill and also some satellites that are stable 
according to Golubev's method are unstable in the sense of Sundman stability. 
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1 INTRODUCTION 

The study of the motion stability of the planet satellites has been usually performed 
by means of the Hill surfaces (Hill, 1878) constructed either for the model of the restricted 
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three-body problem (Proskurin, 1950), or for the problem of Hill (Hagihara, 1952). Since 
Golubev (1967, 1968, 1985) who developed the method it is sometimes referred to as 
c^h method for the general three-body problem, the stability analysis of the motions is 
usually carried out by Golubev's method. This method is based on the famous Sundman 
inequality (Sundman, 1912). 

The search for regions of stable motions in the general three-body problem is divided 
into two tasks: 

- determination of stability regions in the plane of two free parameters — the constant 
of the energy integral and constant of the moment of momentum integral, 

- determination of stable regions in the space of the coordinates used. 

Actually Golubev (1967) carried out the complete solution of the first task by the 
introduction of "the index of Hill stability" s — (?h, where h is the energy constant, and 
c is the moment of momentum constant. Golubev showed that the curve s — Scr is the 
boundary of the stability region in the plane of ch, where the value of Scr is calculated at 
the Eulerian inner libration point L2 with the use of the known equation of fifth power 
as some function of body masses. 

The solution of the second problem is considerably more complex, since it requires the 
construction of the Sundman surfaces in the multidimensional space of the coordinates 
used. For the solution of this problem Golubev (1968) applied the simplified Sundman 
inequality, with the aid of which the solution of problem leads to the construction of Hill 
curves in the plane of two rectangular coordinates x and y. 

Some authors (Marchal, Saari, 1975), (Zare, 1976), (Marchal, Bozis, 1982), (Mar- 
chal, 1990) have elaborated the Golubev's results with the use of the simplified Sundman 
inequality. 

The construction of the regions of possible motions in the space of the selected coor- 
dinates with the use of the exact Sundman inequality is hindered by the large number 
of variables. Thus, in the relative or Jacobian coordinate systems the number of the 
variables is six: three coordinates of one body and three of the other. Therefore, the 
construction of the Sundman surfaces can be carried out in the six- dimensional space of 
these coordinates. 

In the large series of works (Szebehely, Zare, 1977), (Walker, Emslie, Roy, 1980), 
(Donnison, Williams, 1983), (Donnison, 2009), (Li, Fu, Sun, 2010), (Donnison, 2010) 
and other authors the determination of the regions of the possible motions is carried out 
in the six-dimensional space of the Keplerian elements ai,a2, 61,62, ii and ^2 or in the 
four-dimensional space oi, 02, 61 and 62 for the planar problem. 

At the same time, the construction of the Sundman surfaces and, thus, the determina- 
tion of the regions of possible motions and stability regions can be conducted in the space 
of three variables. This substantially facilitates the readabihty of results and eases their 
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application. This possibility is explained by the fact that the exact Sundman inequality 
depends on the coordinates only by means of three quantities — three mutual distances 
between the bodies. 

The article (Lukyanov and Shirmin, 2007) is likely to be the first work that confirmed 
this possibility. In this work the mutual distances between the bodies are the rectangular 
coordinates. The existence of the Hill surfaces's analogue for the general three-body 
problem in the space of the mutual distances is shown here using the exact Sundman 
inequality. The stability regions are determined in the space of the mutual distances and 
have the form of an infinite "tripod". 

Lukyanov (2011) used another choice of three coordinates. The coordinate system 
is determined by the accompanying triangle of mutual positions of bodies. Namely, the 
origin of the coordinate system coincides with one of the bodies, the axis x is directed 
towards the second body, the axis y is perpendicular to the axis x and lies in the plane 
of the triangle. In this coordinate system the position of all bodies is determined by 
three coordinates x, y, R, where x and y are the coordinates of the third body, and R is 
the distance between the first and second bodies. The system of coordinates xyR is to 
a certain degree similar to the rotating coordinate system in the restricted three-body 
problem following the motion of the basic bodies. 

Then, in the space of the coordinates xyR the Sundman surfaces are constructed, the 
singular points of surfaces (coinciding with the Euler and Lagrange libration points) are 
located, the regions of the possible motions and Sundman stability regions are determined. 
The stability region of any body relative to the other body in the space xyR has the form, 
similar to an infinite " spindle" . No restriction on the masses of bodies or their mutual 
positions is assumed in this case. 

In the present work the method of constructing the regions of the possible motions 
(Sundman lobes) in the general three-body problem is presented. The Sundman stability 
analysis of the natural satellites of the planets is carried out, using the high- precision 
ephemeris of the natural satellites of planets calculated on the web-site " Natural Satellite 
Data Center" (NSDC) (Sternberg Astronomical Institute, Moscow) 
{http : / /www. sai.msu.ru/neb/nss /index. htm). 



2 SUNDMAN SURFACES 

The regions of possible motions of bodies in the general three-body problem are de- 
termined by Sundman inequality 

{U-C)J>B, (1) 
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where the force function U and the barycentric moment of inertia J are determined by 
the expressions: 

_ Gmima Gmamg Gm^mi 

R ' ^ 

-"-12 -^23 -n-si 

^ _ mim2-Ri2 + m2m3i?^3 + m^niiRl^ 
m ' 

C = — /i is the analogue of the Jacobi constant, h is the energy constant, B = is the 
Sundman constant, c is the constant of the integral of area. 

Here: G is the universal gravitational constant, mi,m2,m3 are the masses of bodies, 
m = mi + 1712 + rn^ is the total mass of the system, R12, R23, R31 are the mutual distances 
between the bodies. 

Constants C and B are determined by the initial conditions in the barycentric coor- 
dinate system from the relationships: 



V? vi 

C = -h = U - mi— - m2— - "^3^, 
B = — = -(miri X Vi + m2r2 x V2 + mgrg x V3) 



(4) 



Fj, Vj, (z = 1, 2, 3) are the barycentric state and speed vectors of the bodies. 

The boundary of the region of possible motions can be established if in ([T]) inequality 
is replaced with equality: 

{U-C)J = B. (5) 

This equality determines the equation of the Sundman surface. In the general case 
the mutual distances in (|5]) depend on nine coordinates of three moving bodies, which 
substantially hampers the construction of the Sundman surfaces. The transformation to 
the relative coordinate system makes it possible to reduce the number of coordinates to 
six. However, in this case the construction of the Sundman surfaces should be conducted 
in the six- dimensional space. Furthermore, the number of coordinates can be reduced to 
three if the positions of bodies are determined by the following special coordinates. 

The position of the body M2 relative to the body Mi will be characterized by the 
abscissa R on the axis MiX. we will define the position of the body M3 relative to Mi 
by the rectangular coordinates X and Y in the system of MiXY, which always lies in the 
plane that passes through all three bodies. The positions of the bodies in the coordinate 
system MiXYR are defined by three quantities - coordinates X, Y and R, which allows 
us to construct the Sundman surfaces in the three-dimensional space. 

We will use a dimensionless system of coordinates Mixy making the substitution 

X = Rx, Y = Ry. (6) 

Then the mutual distances between the bodies can be expressed in terms of three 
quantities x, y, R. The Sundman surface equation transforms to the form of the functions 
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of three variables 

S{x,y,R) = {U-C)J = 
r 
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Equation ([7j) allows us to conduct the construction of the Sundman surface in the 
three-dimensional cartesian space of variables xyR. 

The singular points of the Sundman surfaces are determined from the system of three 
algebraic equations 
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From the third equation of this system it is possible to determine the mutual distance 
R between the bodies Mi and M2 in the form of the function of unknowns x, y and 
constant G: 
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Substituting this expression for R to the first two equations of set ([8]), we will obtain 
the system of two equations with two unknowns x and y: 

17121713 17131711 
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The second equation in set flTI]]) can be satisfied in two ways: to set y = or to consider 
as zero the entire coefficient in the brackets with y. The first possibility [y = 0) leads to 
the collinear singular points, the second — to the triangular. 

For y = Owe obtain from set (fTOl) one equation for the determination of the coordinate 
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X of the collinear singular points 



(p{x) = 17111712 H : H -^=^ [m2i7i3[x — 1) + m^mixl — 

fx -1)2 



(11) 



[mim2 + 7n27ns{x — 1)^ + m^7nix^] = 0. 



(x- 1)2 xVx\ 

The derivative (p'{x) is always positive, and the following limits take place: 

lim ip(x) = =Foo, lim ip(x) = ±oo, lim (pix) = ±oo. (12) 

This proves the existence of three real solutions of the equation if{x) = 0, which, in 
their turn, determine three collinear singular points of the family of the Sundman surfaces 
in space xyR. 

T I n Gmim2 Gmamg Gm^mA 

Li = Xi, 0, + + ^ , « = 1,2,3 , 13 



where the coordinates Xi are determined by the numerical solution of equation fill I) . 
But if y 7^ 0, then after simple conversions we obtain two triangular solutions of set 

(HDD 

i?23 = R^i = (14) 
which determine two triangular singular points in the space xyR: 



_ / 1 \/3 mim2 + 77127713 + msmA 
^^■^ - 1^2 ' ^1"' ^ 2C ) ■ 



(15) 



The obtained collinear and triangular singular points correspond to the collinear Euler 
and triangular Lagrange solutions known in the general three-body problem. 

Collinear singular points in the space xyR lie in different planes R = Ri, i.e., Ri ^ Rj, 
and for triangular singular points the following equality is fulfilled: R4 = R5. 

Knowing the coordinates of singular points and constant C, from formula ([7]) the 
values of Sundman constant Bi,B2, -B3 and -B4_5 in all singular points Lj, (7 = 1, 2, . . . , 5) 
are calculated. Constants C and Bi are connected by reciprocal proportion 
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where (r23)j and (7'3i)j are calculated at the singular point Lj. 

The relations (1161) have been established by Golubev (1967) in the c^h method. 

Singular points are the points of bifurcation, in which a qualitative change in the 
shape of Sundman surface occurs. The curves f|T6|) on plane GB are the boundaries of the 
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Figure 1: The stable region (hatched) and the boundaries (Li, L2, L3, L4 antiLs) of different 

topological types of the regions of the possible motion for = = — -. 

9 3 1 

topologically different regions of the possible motion. The curve L2 limits the Sundman 
stable region of the body M3, and this stable region is shaded (Fig{T]). 

The general form of the Sundman surfaces for three bodies with the mass ratio pro- 
portional to 9:3:1 is shown in FigJ2l the section of the Sundman surfaces by the planes 
R = R2 and y = is presented in FiglHJ 

In the general three-body problem, as in the restricted problem, the concept of Hill 
stability is conserved. But, to distinguish it from the restricted problem in the general 
three-body problem, we will call this stability Sundman stability. 

We will call the motion of the body M3 in the general three-body problem stable on 
Sundman if there are such regions of the possible motions, limited by the appropriate 
Sundman surfaces, inside which the body M3 will be always (at any instant of time) 
located at a finite distance from one of the bodies Mi or M2. In other words, the body 
M3 will be an eternal satellite of one of Mi or M2 bodies, while bodies Mi or M2 can be 
at any distance one from the other, including infinite. 
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0.15-, 




Figure 2: General view of Sundman surfaces for — = — = — , C=2 and B = Bo- Surfaces 
B 9 3 1 

are represented in the field of the space limited to planes: x = —3, x = 3, y = —1, y = 3, R = 

0, R = 0.7. 

Criterion of Sundman stability is the inequality 

B > B2, (17) 

where B2 is the value of Sundman constant in the inner Euler libration point L2. The 
fulfillment of this condition guarantees that the body M3 can be: in some "spindly" 
surfaces (see Fig.2, 3) remaining the eternal satellite of a body Mi; or in other "spindly" 
surfaces, remaining the satellite of body M2, or in a remote open oval area, when the 
distance between bodies Mi and M2 remains finite, not exceeding Gmim2/C. This last 
case can be treated as Sundman stability of the relative motion of bodies Mi and M2. 

Thus, for (|T7I) any pair of bodies will have Sundman stability if at the initial instant 
the bodies forming this pair are in one of these regions of stability. The loss of stability 
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Figure 3: The sections of the Sundman surfaces for — - = — - = — - by the plane R = R2 (left), 

9 3 1 
by the plane y = (right). 

(body M3 leaving the "spindly" area) occurs if the value R is close enough to its value 
i?2 at the libration point L2. 

3 SUNDMAN STABILITY 

OF THE PLANET SATELLITES' MOTION 

The analysis of Sundman stability of the motion of all known natural planet satellites 
of the Solar System is investigated with the presented theory. The ephemeris of all planet 
satellites are calculated with the most uptodate theories implemented on the NSDC web- 
site (Natural Satellites Data Center) (Sternberg Astronomical Institute, Moscow, Russia), 
constructed by Emelyanov, Arlo (2008) {http : / /www .sai.msu.ru/neb/nss /index. htm) . 
From these ephemeris constants C, B and B2 were calculated in the barycentric coordinate 
system. Sundman stability was determined from formula (fTTI) . 

For each satellite the construction of the Sundman surface sections by the coordinate 
plane xy was also conducted. The sections are given for the Jovian satellite J6 Himalia 
and J9 Sinope (FigJl]). Himalia's Sundman curve is within the Sundman stability region, 
Sinope's curve is outside. In spite of the location of the Sinope orbit inside the Sundman 
lobe, which corresponds to Sundman constant value B = B2, its energy is sufficiently 
high, and it has a potential capability to leave this lobe (see the dash curve). But this 
does not mean that the satellite will leave the vicinity of the planet without fail. Sundman 
instability means that the Sundman surfaces are open and allow it to leave the vicinity 
of planet. But the Sundman surfaces do not tell if this will occur or not. The same case 
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Figure 4: The sections of the Sundman surfaces by the plane of R = R2 for Jovian satelhtes 
J6 Himalia (left) and J9 Sinope (right). Here: the solid line B2 is the boundary of Sundman 
stability region, the dash line is satellite's Sundman curve, the dotted line is approximate region 

of satellite motions. 

is for the Hill stability. 

Lukyanov (2011) showed the Sundman stability of the Moon motion. The Sundman 
stability results for the rest satellites are given in Tables ([UlSj). The tables also list the 
results of the classical Hill stability. All satellites are located in the order of increasing 
semimajor axes of their orbits around the planet. The relative masses of distant planet 
satellites obtained from satellite photometric observations (Emelyanov, Uralskaya, 2011), 
are taken from the web-site NSDC (jhttp : / / www. sai . msu . ru / neb / nss / index. htm ) . 

The Martian satellites (Phobos and Deimos) have Hill and Sundman stability (TabJT]). 

The main and the inner satellites of Jupiter have Hill stability and Sundman stability 
and are not included in the tables. The distant satellites, which have prograde and 
retrograde orbits, are of special interest. All prograde satellites of Jupiter have Hill and 
Sundman stability (Tabl2]). All retrograde satellites with a > 18.34 ■ 10® km are unstable 
according to Hill and Sundman, independently from their masses. An exception is the 
satellite S/2003 J12 (a = 19 ■ 10® km, i = 145.8°, e = 0.376) with a relatively small mass, 
which has Hill stability and Sundman stability. 

The situation is different for the satellites of Saturn. The main, inner and distant 
prograde satellites of Saturn, which belong to the Gallic group {i = 34°) and Inuit group 
{i = 45°), have Hill stability and Sundman stability (Tab. [3]). The retrograde satellites 
with a < 18.6 ■ 10® km have Hill and Sundman stability, with a > 18.6 ■ 10® km have Hill 
stability, but Sundman instability. Furthermore, Sundman unstable is also the satellite S 
LI Greip with the semimajor axis a = 18.1 ■ 10® km. 

The main and inner satellites of Uranus have the Hill stability and Sundman stability. 
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The stability results coincide for all distant satellites, except for the most distant satellite 
U XXIV Ferdinand (a = 20.9 ■ 10^ km), which has Sundman instability (TabJl]). 

Triton and the Neptune inner satellites have the Hill and Sundman stability. Two 
distant Neptune satellites have Sundman instability: N X Psamathe (a = 46 ■ 10® km) 
and N XIII Neso (a = 48 ■ 10® km). The rest of Neptune satellites have Hill stability and 
Sundman stability (Tab ED- 

The comparison of the results of the Sundman stability and Hill stability shows that 
Hill stability always follows from the Sundman stability, but the reverse assertion is not 
correct. It is caused by the fact that, in contrast to Hill's model, in the Sundman model 
the satellite masses are not zero, but are finite. Therefore, each satellite of any planet has 
an individual value of the Sundman constant B2, while in Hill's model all satellites of any 
planet have the same value of the Hill constant C2. 

The comparison of the obtained results with Golubev c^h method is carried out in two 
directions: 

— comparison of the stability criteria used, 

— comparison of the obtained regions of the possible motions. 

Analytical forms of stability criterion in our work B > B2 and in Golubev's method 
c^h < {c'^h)cr, are the same. But the calculation of the constants in the left and right sides 
of the inequalities is carried out using different formulas. This leads to some differences in 
numerical results. The comparison with the results of the work (Walker et al, 1980) for the 
satellites J1-J13 shows that the Sundman stability or instability of the these satellites, 
obtained in our work, agrees with the results of Walker et al. (1980) for all satellites, 
except for four satellites of Jupiter with retrograde motion, J VIII, J IX, J XI and J XII. 
For these satellites we obtained instability, while in the work cited these satellites were 
indicated as being stable. This is likely to be due to the approximation of the three-body 
problem by two problem of two bodies and also by the neglect of orbit inclinations. 

We conducted the construction of regions of possible motions in the three-dimensional 
space of xyR, while in all works of other authors the value of R is excluded from the 
examination, and the constructions of regions of possible motions are conducted in the 
xy-plane. For this reason in the c^h method it is not possible to get a number of important 
results. For example, it cannot be obtained that the loss of stability (withdrawal of the 
body M3 from the stability region) can occur only for a certain distance between the 
bodies Mi and M2. Generally, Sundman curves in the plane R = const with a change in 
R can sharply and qualitatively differ from Hill's curves, as shown by Lukyanov (2011). 
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4 DISCUSSION 



The famous Sundman inequality in the general three-body problem takes the form 

iU-C)J-B>^. (18) 

For the material motions of bodies, i.e., with the fulfillment of conditions > 0, it 
determines the regions of possible motions satisfying the inequality 

{U-C)J>B. (19) 

The boundaries of the region of possible motions are determined by the equation 

{U-C)J = B, (20) 

which we call the equation of the Sundman surface, while the stability in Hill's sense for 
the three-body problem — Sundman stability. By analogy with the surfaces of the zero 
speed in the restricted three-body problem, we may call the Sundman surfaces in the 
general three-body problem the surfaces of zero rate of change of the barycentric moment 
of inertia of bodies ( J = 0). 

The determination of the Sundman stability and the construction of the Sundman 
curves in the plane of parameters C and B (see Fig. 1) is completely solved by Golube^|3 
(1967) in his c'^h method (in our designations = 2B,h = —C). Now this method is 
called Golubev's method. 

Golubev's method determines not the surfaces, but the Sundman curves located in 
the plane of the triangle formed by the mutual distances between the bodies. The mutual 
distances between the bodies R13 and R23 are substituted by the relative values R13/R 
and R23/R, and the value of i? = R12 is generally excluded from examination. 

Equation of "current" Sundman curve in Golubev's method has the form of the hy- 
perbola CB = const. If in this case the constants C and B are expressed in term of any 
other variables, then, in its turn, the task of construction of the Sundman curves in the 
space of these variables arises. Thus, in the large series of works of (Szebehely and Zare, 
1976), (Walker, 1983), (Donnison, 2010) and many other authors the task of construct- 
ing Hill-Sundman curves and determination of stability regions in the general three-body 
problem is solved by Golubev's method in the space of six quantities: semimajor axes 
ai, 02, eccentricities ei, 62 and inclinations ii,i2, for calculation of the constants C and B 
the approximation of three-body problem by two problems of two bodies is used. This 
introduces a certain error to the solution of problem. Besides, the value of R remains 
unknown. 



^in the English-language literature the surname Golubev is frequently written incorrectly. 
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For the representation of the Sundman curves on the plane xy, Golubev (1968) con- 
sidered another method. He used the simphfied Sundman inequahty instead of the exact 
inequahty f lT9|) 

f/V > BC, (21) 

which is the consequence of inequahty (I19p and is obtained after the multiphcation of 
inequahty ( IT^ by f/, taking into account inequahties C > and U > C. Inequahty ([21]) 
does not reflect the entire diversity of the Sundman surfaces. 

Like the c^h criterion (obtained from the condition of the positivity of the discriminant 
of the quadratic trinomial for R from the left side of the Sundman inequality), simplified 
inequality fl^T]) does not contain the mutual distance R12 = R- Therefore, by means of 
inequality ( l2T|) . it is possible to construct not the surfaces, but the Sundman curves in 
the plane of relative coordinates xy. The construction of these curves was subsequently 
conducted in the works of (Marchal, Saari, 1975), (Marchal, Bozis, 1982) and other 
authors. 

Thus, the task of constructing the Sundman surfaces in the space of the coordi- 
nates used remained incomplete before the publication (Lukyanov, Shirmin, 2007) and 
(Lukyanov, 2011) appeared. Lukyanov, Shirmin (2007) used the mutual distances be- 
tween the bodies as the coordinates. This made possible to construct exact Sundman 
surfaces in the three-dimensional space of mutual distances. Lukyanov (2011) used the 
more convenient rectangular coordinate system xyR, determined by the accompanying 
triangle of mutual positions of three bodies. 

In these works the exact Sundman inequality (fT9l) is used and, therefore, the value of 
R is not excluded from the examination. In this case no simplifications or assumptions 
are applied. The construction of the Sundman surfaces is implemented in the three- 
dimensional space of the coordinates used with the determination of the singular points 
of surfaces, regions of the possible motion and Sundman stability regions. 

Regions of the possible motion constructed by means of the exact Sundman inequalities 
differ from analogous regions defined according to the simplified Sundman inequality, both 
quantitatively and qualitatively. 

The stability regions determined by the simplified Sundman inequality (l2Ti) have larger 
sizes than those calculated by exact inequality iHM . Therefore, the stability obtained by 
means of (ETi) can turn to instability, when using exact inequality (IT^ . 

It is easy to derive by means of the exact Sundman surfaces that the loss of Sundman 
stability for the body M3 can occur only when a certain distance R between the bodies 
Ml and M2 takes place, so that the "passage" through the neighborhood of the singular 
point L2 is open. It is caused by the fact that the singular points of the Sundman surfaces 
are determined by three coordinates Li{xi, yi, Ri) and in the space xyR they lie, generally 
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speaking, in different planes. This result cannot be established with the aid of inequality 
(^I^, since it does not depend on R. 

The construction of exact Sundman surfaces allows us to define the regions of possible 
motions for any of the three bodies and for any values of C and B. Using the Sundman 
surfaces yields, for example, that with the fulfillment of the stability criterion the body 
M3 (it can be any body) for any time — 00 < t < 00 will be located at a finite distance 
from one of the bodies Mi or M2 or at a large distance from these bodies. Qualitatively, 
the analogous result is known for the Hill surfaces in the restricted three-body problem 
as well. If the body M3 is located, for example, in the stability region near Mi, then the 
Sundman surfaces admit the possibility of retreating of the body M2 to any large distance 
from the pair Mi, M3. For the Hill surfaces this situation is not possible. 

By means of the Sundman surfaces it is possible to establish the stability of only 
one pair of bodies, and the third body will be in this case unstable in the Sundman 
sense. Sundman surfaces do not establish the simultaneous stability of three bodies, i.e., 
guaranteed location of all bodies in a certain finite region of the space (Lagrange stability), 
although these surfaces do not exclude this case. Sundman instability does not mean that 
a body will necessarily leave the neighborhood of another body. The Sundman surfaces 
do not allow us to determine if this retreat will actually occur. This result is analogous to 
that of Hill stability. The determination of Sundman stability of the planet satellites of 
the Solar system conducted in this study shows the effectiveness of the use of Sundman 
surfaces in the coordinate form. 

We believe that our results represent a certain interest for celestial mechanics and for 
astronomy clS db whole. 
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1: Martian satellites. Here a is the semimajor axis of the satellite orbit, i is the inclinati 
e is the eccentricity, ra/Mp is the ratio of the satellite mass to the planet mass. 



Satellite 


a 


e 


i 


m/Mp 


Stability 




{km) 




(deg) 


10-8 


Hill 


Sundman 


Ml Phobos 


9380 


0.0151 


1.1 


1.6723 


yes 


yes 


M2 Deimos 


23460 


0.0002 


0.9 - 2.7 


0.2288 


yes 


yes 



Table 2: The irregular Jovian satellites (the notation in Tabled]). 



Satellite 


d 


i 




g 


m/Mp 


Stability 




(lO^km) 

\ \J III 1 1 V 1 


idea') 






10-9 


Hill 


SmITI mTTIPI Tl 
kJ LliiU-liiCXii 


I 


2 






4 






7 


yv V 111 J- lidiiioiu 


7 ^n? 




n 


942 




yes 


yes 


XIII Leda 


n 1 65 


27 46 


n 


164 


5 76 


yes 


yes 


\/T l~ri m PI 1 1 PI 

V -L -L XlXliCXllCL 


11.461 


27 50 


n 


.162 


22101 8 


yes 


yes 




1 1 71 7 




n 


1 1 9 
.11^ 


OO 1 . tj 


yes 


yes 


VTT F.larp 


1 1 741 
11. 1^1 


9fi fi'^ 


n 


91 7 


4^7R 9 


yes 


yes 


XTVT riarnn 


1 6 989 


^1 4 


n 


430 


3394 


yes 


yes 


S/9003 13 


1 8 340 


143 7 

1t:0. 1 


n 


941 


1 963 


no 


no 


S/9003 11 9 


1 9 009 


14^ 8 


n 


376 


0631 


yes 


yes 


XXXIV Eiinorie 


1 9 302 


145 8 


n 


1 44 


9447 


no 


no 


S/2003 J18 


20.700 


146.5 


0, 


.119 


0.2920 


no 


no 


XXXV Orthosie 


20.721 


145.9 


0, 


.281 


0.3315 


no 


no 


XXXIII Euanthe 


20.799 


148.9 


0, 


.232 


0.4341 


no 


no 


XXIX Thyone 


20.940 


148.5 


0, 


.229 


0.6946 


no 


no 


S/2003 J16 


21.000 


148.6 


0, 


.270 


0.1342 


no 


no 


XL Mneme 


21.069 


148.6 


0, 


.227 


0.3315 


no 


no 


XXII Harpalyke 


21.105 


148.6 


0, 


.226 


0.8367 


no 


no 


XXX Hermippe 


21.131 


150.7 


0, 


.210 


1.4919 


no 


no 


XXVII Praxidike 


21.147 


149.0 


0, 


.230 


2.8495 


no 


no 


XLII Thelxinoe 


21.162 


151.4 


0, 


.221 


0.3473 


no 


no 


XXIV locaste 


21.269 


149.4 


0, 


.216 


1.3971 


no 


no 


XII Ananke 


21.276 


148.9 


0, 


.244 


157.9 


no 


no 


S/2003 J15 


22.000 


140.8 


0, 


.110 


0.1342 


no 


no 


S/2003 J4 


23.258 


144.9 


0, 


.204 


0.0947 


no 


no 


L Herse 


22.000 


163.7 


0, 


.190 


0.2526 


no 


no 



Continued on the next page 
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Tables 2 continued. 



1 




2 


3 






4 


5 


6 


7 


S/2003 J9 


22, 


.442 


164. 


,5 


0, 


.269 


0.0947 


no 


no 


S/2003 J19 


22, 


.800 


162. 


,9 


0, 


.334 


0.1263 


no 


no 


XLlll Arche 


22, 


.931 


165. 


,0 


0, 


.259 


0.2842 


no 


no 


XXXVIII Pasithee 


23, 


.096 


165. 


,1 


0, 


.267 


0.1658 


no 


no 


XXI Cnaldene 


23, 


.179 


165. 


,2 


0, 


.251 


0.7499 


no 


no 


XXXVII Kale 


23, 


.217 


165. 


,0 


0, 


.260 


0.2447 


no 


no 


XXVI Isonoe 


23 


.217 


165. 


,2 


0, 


.246 


0.6157 


no 


no 


XXXI Aitne 


23, 


.231 


165. 


,1 


0, 


.264 


0.4026 


no 


no 


XXV Erinome 


23, 


.279 


164. 


,9 


0, 


.266 


0.3789 


no 


no 


XX Taygete 


23, 


.360 


165. 


,2 


0, 


.252 


1.1445 


no 


no 


XI Carme 


23, 


.404 


164. 


,9 


0, 


.253 


694.6 


no 


no 


XXIII Kalyke 


23, 


.583 


165. 


,2 


0, 


.245 


1.5471 


no 


no 


XLVII Eukelade 


23, 


.661 


165. 


,5 


0, 


.272 


0.7104 


no 


no 


XLIV Kallichore 


24, 


.043 


165. 


,5 


0, 


.264 


0.2289 


no 


no 


S/2003 J5 


24, 


.084 


165. 


,0 


0, 


.210 


0.9788 


no 


no 


S/2003 JIO 


24, 


.250 


164. 


,1 


0, 


.214 


0.0947 


no 


no 


XLV Helike 


21, 


.263 


154. 


,8 


0, 


.156 


0.7183 


no 


no 


XXXII Eurydome 


22, 


.865 


150. 


,3 


0, 


.276 


0.4262 


no 


no 


XXVIII Autonoe 


23 


.039 


152. 


,9 


0, 


.334 


0.7814 


no 


no 


XXXVI Sponde 


23, 


.487 


151. 


,0 


0, 


.312 


0.2763 


no 


no 


VIII Pasiphae 


23, 


.624 


151. 


,4 


0, 


.409 


1578.7 


no 


no 


XIX Megaclite 


23 


.806 


152. 


,8 


0, 


.421 


2.1312 


no 


no 


IX Sinope 


23, 


.939 


158. 


,1 


0, 


.250 


394.7 


no 


no 


XXXIX Hegemone 


23, 


.947 


155. 


,2 


0, 


.328 


0.3394 


no 


no 


XLl Aoede 


23, 


.981 


158. 


,3 


0, 


.432 


0.6473 


no 


no 


S/2003 J23 


24, 


.055 


149. 


,2 


0, 


.309 


0.0947 


no 


no 


XVII Callirrhoe 


24, 


.102 


147. 


,1 


0, 


.283 


5.3044 


no 


no 


XLVIII Cyllene 


24, 


.349 


149. 


,3 


0, 


.319 


0.2368 


no 


no 


XLIX Kore 


24, 


.543 


145. 


,0 


0, 


.325 


0.3947 


no 


no 


S/2003 J2 


28, 


.570 


151. 


,8 


0, 


.380 


0.1500 


no 


no 
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Table 3: The irregular Saturnian satellites (the notation in Tabled]). 



kJClljt'lllljt' 




Li 


b 






m 1 A/f r. 


Stability 




(10 




{deg) 






iU 


Hill 


Sundman 


XXIV Kiviuq 


11, 


.111 


45.71 


0. 


,334 


0.8629 


yes 


yes 


XXII Ijiraq 


11, 


.124 


46 .U 


0, 


.316 


0.3248 


yes 


yes 


IX Phoebe 


12, 


.944 


174.8 


0, 


,164 


1458.957 


yes 


yes 


XX Paaliaq 


15, 


.200 


45.13 


0, 


.364 


2.2728 


yes 


yes 


XXVII Skathi 


15, 


.541 


152.6 


0, 


,270 


0.0588 


yes 


yes 


XXVI Albiorix 


16, 


,182 


33.98 


0, 


,478 


4.3629 


yes 


yes 


S/2007 S2 


16, 


.560 


176.7 


0. 


,218 


0.0248 


yes 


yes 


XXXVII Bebhionn 


17, 


.119 


35.01 


0, 


.469 


0.0261 


yes 


yes 


XXVIII Erriapus 


17, 


.343 


34.62 


0, 


.474 


0.2294 


yes 


yes 


XXIX Siarnaq 


17, 


.531 


45.56 


0, 


.295 


24.1988 


yes 


yes 


XLVII Skoll 


17, 


.665 


161.2 


0, 


.464 


0.0237 


yes 


yes 


LII Tarqeq 


17, 


.920 


49.86 


0, 


.107 


0.0385 


yes 


yes 


XXI Tarvos 


17, 


.983 


33.82 


0, 


,531 


0.5455 


yes 


yes 


LI Greip 


18, 


.105 


172.7 


0, 


.374 


0.0158 


yes 


no 


XLIV Hirrokkin 


18, 


.437 


151.4 


0, 


,333 


0.0965 


yes 


yes 


S/2004 S13 


18, 


.450 


167.4 


0, 


,273 


0.0148 


yes 


yes 


S/2004 Sir 


18, 


.600 


166.6 


0, 


,259 


0.0082 


yes 


yes 


L Jarnsaxa 


18, 


.600 


162.9 


0, 


.192 


0.0116 


yes 


no 


XXV Mundilfari 


18, 


,685 


167.3 


0, 


,210 


0.0464 


yes 


no 


S/2006 SI 


18, 


.981 


154.2 


0, 


,130 


0.0192 


yes 


no 


XXXI Narvi 


19, 


.007 


145.8 


0, 


.431 


0.0340 


yes 


no 


XXXVIII Bergelmir 


19, 


.338 


158.5 


0, 


.142 


0.0248 


yes 


no 


XXIII Suttungr 


19, 


.459 


175.8 


0, 


.114 


0.0422 


yes 


no 


S/2004 S12 


19, 


.650 


164.0 


0. 


,401 


0.0142 


yes 


no 


S/2004 S07 


19, 


.800 


165.1 


0. 


,580 


0.0200 


yes 


no 


XLIII Hati 


19, 


.856 


165.8 





.372 


0.0185 


yes 


no 


XXXIX Bestla 


20 


.129 


145.2 


0. 


,521 


0.0432 


yes 


no 


XL Farbauti 


20, 


.390 


156.4 


0. 


,206 


0.0113 


yes 


no 


XXX Thrymr 


20, 


.474 


176.0 


0, 


,470 


0.8278 


yes 


no 


S/2007 S3 


20, 


.518 


177.2 


0, 


,130 


0.0119 


yes 


no 


XXXVI Aegir 


20, 


.735 


166.7 


0, 


,252 


0.0214 


yes 


no 


S/2006 S3 


21, 


.132 


150.8 


0, 


,471 


0.0100 


yes 


no 


XLV Karl 


22, 


.118 


156.3 





.478 


0.0409 


yes 


no 


XLI Fenrir 


22, 


.453 


164.9 


0, 


,136 


0.0095 


yes 


no 


XLVIII Suit 


22 


.707 


177.5 





.451 


0.0127 


yes 


no 


XIX Ymir 


23, 


.040 


173.1 





.335 


1.3878 


yes 


no 


XLVI Loge 


23, 


.065 


167.9 


0, 


.187 


0.0232 


yes 


no 


XLII Fornjot 


25, 


.108 


170.4 


0. 


,206 


0.0211 


yes 


no 
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Table 4: The irregular Uranus' satellites (the notation in Tabled]). 



Satellite 


a 


e 


i 


m/Mp 


Stability 




(lO^km) 




(deg) 


10-9 


Hill 


Sundman 


XXII Francisco 


4.2760 


0.1425 


147.613 


0.0658 


yes 


yes 


XVI Caliban 


7.1689 


0.0823 


139.681 


8.1305 


yes 


yes 


XX Stephano 


7.9424 


0.1459 


141.538 


0.3494 


yes 


yes 


XXI Trinculo 


8.5040 


0.2078 


166.332 


0.0593 


yes 


yes 


XVII Sycorax 


12.2136 


0.5094 


152.669 


46.6790 


yes 


yes 


XXIII Margaret 


14.3450 


0.7827 


50.651 


0.0609 


yes 


yes 


XVIII Prospero 


16.1135 


0.3274 


146.340 


1.1306 


yes 


yes 


XIX Setebos 


18.2052 


0.4943 


148.828 


1.4240 


yes 


yes 


XXIV Ferdinand 


20.9010 


0.4262 


167.278 


0.0874 


yes 


no 



Table 5: The irregular Neptune's satellites (the notation in Table [TJ. 



Satellite 


a 


e 


i 


m/Mp 


Stability 




(lO^fcm) 




{deg) 


10-9 


Hill 


Sundman 


II Nereid 


5.5134 


0.7512 


7.232 


301.38 


yes 


yes 


IX Halimede 


15.728 


0.5711 


134.101 


3.0835 


yes 


yes 


XI Sao 


22.422 


0.2931 


48.511 


0.6445 


yes. 


yes 


XII Laomedeia 


23.571 


0.4237 


34.741 


0.5606 


yes 


yes 


X Psamathe 


46.695 


0.4499 


137.391 


0.9244 


yes 


no 


XIII Neso 


48.387 


0.4945 


132.585 


1.3423 


yes 


no 
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